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TECHNICAL NOTE 4408 


THE THEORY OF DIFFUSION IN STRAINED SYSTEMS 
By Louis A. Girifalco and Hubert H. Grimes 


SUMMARY 

Because the current theory of solid-state diffusion is limited to 
unstrained crystals and cannot be applied readily to strained systems. 
Fief a first and second lavs vere generalized to include the effects of 
strain on the diffusion rates. The nonhomogeneity introduced into the 
atomic jump frequency by strain vas found to contribute strain- dependent 
terms to the diffusion equations in addition to the terms containing the 
concentration gradient. 

From a consideration of the effect of strain on the free energy of 
activation, it can be shown that for simple strains, such as those re- 
sulting from compression, tension, shear, and hydrostatic pressure, the 
diffusion coefficient is an exponential function of the lattice parameter. 
An examination of the available experimental data for the variation of 
diffusion coefficients with pressure confirms this theoretical prediction. 

The theory presented herein states that the magnitude of the varia- 
tion of the diffusion coefficient with pressure depends on the interatomic 
forces as the diffusing atom moves from its equilibrium position to the 
activated position. On the basis of this theory, a parameter depending 
upon the interatomic forces can be computed from the experimental data. 

In all cases investigated, the magnitudes of this parameter were in agree- 
ment with the known characteristics of the interatomic potential-energy 
functions of the systems. 

The effect of plastic flow on the diffusion rate was also studied 
by considering the rate at which vacancies are produced by dislocation 
motion and the rate at which vacancies condense at inhomogeneities In 
the crystal. The resulting equations predict that for a vacancy mecha- 
nism the diffusion coefficient varies linearly with the strain rate. 

This conclusion is in agreement with experiment. 
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INTRODUCTION 

The theory of diffusion in solids has been the subject of a great 
deal of investigation in recent years and satisfactory theoretical models 
have been constructed that adequately describe the basic diffusion proc- 
esses in many simple solids. Present theories, however, are limited to 
unstrained crystals and are not strictly applicable to strained systems. 
Since the diffusion rate is determined by the energy of interaction be- 
tween the diffusing atom and the crystal lattice, and since this energy 
depends on the interatomic distances, it Is to be expected that the dif- 
fusion coefficients will' be altered by a strain superimposed on the crys- 
tal. Some experimental evidence is available that indicates that elastic 
strain can increase the self-diffusion coefficient by as much as a factor 
of 2 (ref. l) and that plastic strain can increase the self- diffusion 
coefficient by an order of magnitude (refs. 2 and 3). Also, it is well 
known that hydrostatic pressure decreases the diffusion coefficient; in 
fact, a pressure of 7500 atmospheres is sufficient to lower the self- 
diffusion coefficient of sodium by an order of magnitude at 90° C 
(ref. 4) . 

If the crystal is strained in a nonhomogeneous manner, another fac- 
tor becomes operative in addition to those that control the change In the 
diffusion coefficient. According to the theory of irreversible processes 
(ref. 5), every thermodynamic flux is proportional to every thermodynamic 
force so that the diffusion flux is proportional not only to the concen- 
tration gradients, but Is also proportional to the strain gradients in 
the crystal. Thus, not only is the magnitude of the diffusion coefficient 
changed by a generalized strain, but the basic character of the diffusion 
equations is also changed. 

Because of the important role played by diffusion processes In many 
solid-state phenomena such as oxidation, the annealing of radiation damage, 
creep, and rupture, and in view of the wide variety of applications in 
which materials are under strain, a thorough understanding of the effects 
of strain on diffusion is highly desirable. 

In this report the theory of diffusion in strained systems Is developed 
from the point of view of molecular kinetics, the fundamental physical fac- 
tors involved are discussed, and the resulting theory is compared with ex- 
isting experimental data. 


BASIC EQUATIONS OF DIFFUSION THEORY 

Generalizations of Fick's first and second laws are obtained for the 
diffusion of a single species in a crystalline lattice by a modification 
of the method of conditional transition probabilities (ref. 6). The dif- 
fusion equations are expressed in terms of atomic jump frequencies without 
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the usual condition that the jump frequency is independent of position. 

In this form, the equations can easily be applied to a strained lattice. 

Consider a volume element d? in the crystal centered about the 
point defined by the position vector ?. The number of particles of the 
diffusing species contained in d ? is given by N(?,t)d? where E(?,t) 
is the concentration of diffusing species at the position r and time t. 
In general, the nu m ber of particles in the volume element dr is not 
constant because particles are continually jumping out of d? while other 
particles are jumping into d? from adjacent portions of the crystal. 

The rate at which N(?, t) changes as a result of these two processes can 
be calculated as follows: If A(r,r*,t) dr'dt is the conditional proba- 

bility that an atom in the volume element dr jumps to the volume element 
dr* during time dt, then, the number of particles that jump from dr to 
d?‘ in time dt is given by 

N(?,t) A(r,?’ ,t)dr dr* dt (l) 

and the total number of particles that jump out of d? during time dt 
is obtained by integrating over- all dr* ; that is, 

N (r,t)dr’ dt = N(r, t)dr dt / A(r,r‘ ,t)d?' (2) 

JP 

where N (r,t) is the rate at which particles leave the volume element 
dr. 


Similarly, the number of particles that jump from dr’ to dr dur- 
ing time dt is given by 

K(r,t) A(r’ T ,r,t)dr dr 1 dt (3) 

and the total n umb er of particles that jump into d? from other parts of 
the crystal is given by 

H + (r,t)dr dt = dr dt H(r' , ,t) A(r‘ ,r,t)dr 1 (4) 

where K + (r,t) is the rate at which particles enter the element dr. The 

net rate of increase of N(r,t) is obtained by subtracting equation (2) 
from equation (4) : 


5t = 



N(r',t) A(r* ,r, t)!? 1 - N(?,t) 



A(?,?*,t)d? ! 


( 5 ) 
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At this point it is convenient to express- equation (5) in terms of 
the Jump distance (r*-r) by performing a transformation of variables so 
that 

X = r ' - r (6) 
In terms of the Jump vector X, equation (5) becomes 


St = 



K(r + X,t) A(r + X,r,t)dX - N(r,t) 



A(r,r + X,t)dX 


(7) 


The function A(r" + }£,r,t) is the probability frequency that a par- 

tide at r + X will jump to r, and if the vector X would always 
terminate at a point that is capable of accommodating the diffusing par- 
ticle, A would equal the atomic jump frequency. There are cases, how- 
ever, such as in diffusion by the vacancy mechanism, in which the terminus 
of cannot always accommodate the migrating particle. The jump fre- 
quency must then be multiplied by the probability that a^vacant site 
exists at the end of the atomic jump. Therefore, if T(r,X) is the atomic 
jump frequency with jump vector >£ and n(r",t) is the probability that 
a site at ■? is vacant (i.e., it can accept the jumping particle), then. 


A(r + X,r,t) - T(x + X, -X) n(r,t) 

A(r,r + %t) = r(r,+X) n(r" + 5f,t) 

Substituting equation (8) into equation (7) and arbitrarily replacing X 
with the equally valid negative argument -it in the first integral of 
equation (7) gives 



St = 



k 


N(r“- X,t) r(p - X,X)dX - H(r,t) f T(r,+X) n(r + X,t)df 

Jx 


(9) 


In crystalline solids, T is zero for all jump vectors except for a 
small number of X's, which may be denoted by T may then be expressed 

as a delta function 

r(r,x) = ^ - ^i> (10) 




vhere a is the total number of possible jumps a particle can make out 
of its position at r. The set of possible jump vectors and the 


r 
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value of a are determined by the crystal structure of the lattice. 
Substituting equation (10) into equation (9) and making use of the proper- 
ties of the delta function yields 



x) s(x - x ± )ax 





a 


Y ri(?) 


.-f —b . .*+ “► . -► 

6(X - X ± ) n(r + X,t)dX 


= n(r,t) N(r - X-pt) r ± (r - 5^) 
i=l 


-N(r,t) r ± (r) n(r + jt^t) (ll) 

It will now be assumed that the functions Nr and n can be ex- 
pressed as a Taylor expansion in powers of X^ about the point r, re- 
taining only the first three terms of the expansion. Therefore, 

N(r - %±,t) Ij^r* - ^_) = N(ivt) r ± (r) + % • VNI^ + ^ (XjL • V) 2 NI^ 

( 12 ) 

and 

n(r + ^tpt) = n(r,t) + ^ • Vn + ^ (lt^ • V) 2 n (13) 


W 
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The derivatives are evaluated at the point Equation (ll) can 

now be written as 

= n(r,t) N(r,t) ^ T ± (r) + n(r,t) ^ • VNr ± + 


\ n(r,t) ^ (X ± • V) 2 K r - 


i=l 


-At) *,*) V Ii(r) - K(r,t) r. Ii(r) ? ± . Vn - 


or 


a 

f H(r,t) ^~ , ^ r 1 (r)(t jL • V ) 2 n 
1=1 


(14=) 


ff - 2 |-(\l • ^(nKTi) + | [n(^i • V) 2 Nr ± - HT^ • V) 2 n]J 

(15) 

Equation (15) is a generalization of Fids’ s second law and is valid 
for any system regardless of the nature of or of the coordinate 

system chosen. Also, if the strain is homogeneous, VP^ = 0. With these 
restrictions, equation (15) reduces to 

|| = \ n(r,t)7> l( *i • V) 2 If - \ N(r,t) ^ • V ) 2 n 

1=1 i==1 (16) 

The position vector x in equation (16) is referred to an arbitrary 
coordinate system. It Is always possible to find a transformation of 
coordinates so that the position vector is referred to the principle axes 

of diffusion in which mixed derivatives, such as d^/cbc Sy, vanish. 
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In this case, 


a 



V (\ 2 2 ^2 n 2 52 n \ 

r A Xix + x iy + x iz ) 


N(r,t) \ ' r L 2 d 2 n > 2 d 2 n ,2 d 2 n\ 

2 2 _, H lx 3x2 - X iy 3y 2 + ^ Sz2 ; 


(17) 


where X^ x , Xjy, and X^ z are the components of X^. The coordinates 
are now referred to the principle axes of diffusion. 


Since Pick’s first and second laws are connected by the equation of 
continuity. 


Sn 

St “ - 


7 



(18) 


where J N is the flux of species N, Fick’s first law corresponding to 
the second law given by equation (17) is 


a a 

■r -n(r,t) V ^ r ,2 dH Sf(r,t) \ ^ r , 2 dn 

J x = ~ A 2 / , r i X ix ^ + - g-/ , ^x ^ 

i=l i=l 

J y - i + I 

i=l i=l 


T -n(?,t) \ ^ F 1 2 SlT NCrit) \ ,2 dn 

J z = "2 Zu 1 iz ^ 2~/ /Az ^ 

i=l i=l 

where J x , Jy, and J z are the components of ?jj. 


(19) 


The physical significance of the theoretical development up to this 
point can be seen most clearly by a consideration of equation (14) . In 
the usual expression for Fick’s second law in unstrained systems, al 1 
terms except those involving the second derivative of N are zero. The 
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other terms appear in equation (14), however, because it was not assumed 
in the present development that F ^ and n are independent of position. 

The third term does not sum^to zero since, in general, it is not composed 
of terms antisymmetric in X^. Furthermore, since and n are func- 

tions of position, they must be retained in the differential operators, 
and terms dependent on the first and second derivatives of and n 

appear in equation (14) . 

In the case of a uniform homogeneous strain, equation (14) reduces 
to equation (17). If the gradient of n is zero, equation (17) becomes 
equivalent to Fick’s second law. In this simple case, it is necessary 
only to calculate the effect of the homogeneous strain on and n in 

order to specify the effect of strain on diffusion. The effect of strain 
on the jump frequency can be analyzed in terms of rate theory. This analy- 
sis is presented in the following section. The effect of strain on n is 
dependent upon the diffusion mechanism. For interstitial diffusion, n = 1 
provided that the concentration of interstitials is low and that strains 
do not affect this value . For diffusion by a vacancy mechanism, however, 
n is the vacancy concentration and will vary with Btrain. An analysis 
of this variation is presented in the section DEPENDENCE OF VACANCY CON- 
CENTRATION ON STRAIN. 


DEPENDENCE OF JUMP FREQUENCY ON STRAIN 


According to the statistical theory of rate processes, the jump fre- 
quency is determined by the ratio of two configurational partition func- 
tions, one referring to the activated state and the other referring to 
the normal state. In analyzing the effect of strain on the jump frequency, 
the formulation of the rate process theory in solids given in reference 7 
Is used In which the jump frequency Is given by 


r = 



I 

I 


exp (-cp/kT) da 


exp (-cp/kT) dv 


( 20 ) 


where k is Boltzmann’s constant, T is the temperature, and <p is the 
potential energy of the system as a function of all the coordinates of 
all the atoms in the crystal. The integral in the numerator of equation 
(20) Is evaluated over a hyper surface a in the configuration space so 
that the surface passes through the point corresponding to the diffusing 
atom at its activated position with all other atoms at their equilibrium 
position. The hypersurface Is also required to be perpendicular to 
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contours of constant potential energy in the configuration space. The 
hypersurface defined in this manner divides the configuration space into 
two symmetric parts. The integral in the denominator is evaluated over 
the configuration volume v of one of these symmetric parts. 

Equation (20) was derived for the case of an unstrained crystal. 
However, it is applicable to strained crystals if the potential energy 
q> is taken to be a function of the six strain components e a p as well 

as the atomic .coordinates q* . A similar procedure has been used in ref- 
erence 8 in an analysis of the statistical mechanics of crystal lattices. 
Thus, the potential energy in equation (20) is given by 

<P - <P(< e ap ) (21) 

where q^ represents the set of al 1 atomic coordinates and e a g repre- 
sents the set of six independent strain components. 

The potential* tp can be expanded in a Taylor series about the nor- 
mal lattice configuration in the unstrained state to give 





where q^ represents the set of atomic coordinates when all the atoms 

are in their mean positions for the normal state of the crystal, and the 
double zero subscript indicates that the derivatives are evaluated at the 

point (qj,0). 


W 
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If the potential energy is expanded about the point (q?, 0), where 
+ « 
qT represents the mean atomic coordinates when the system is in the 

activated state, then, 


<p = <p(q.*>o) 




8qj6q k + 


\ 

a,p 

r,8 


6 ap 6 yS 


+ 


2 

J*a,p 


8<1 ' ]6aP (seapO*,0 


(23) 


The subscript t , 0 indicates that the derivatives are evaluated at the 
point (qj,0). 


Since the point (qj,0) corresponds to an unstrained crystal with all 

atoms in a mean position so that the crystal Is at the bottom of a poten- 
tial well, the first derivatives in equation (22) vanish: 



(24) 


The point (q^,0) corresponds to a saddle point with respect to the 
<J 

atomic coordinates q., bo that the derivative of <p with respect to 

J + 

qj also vanishes at this point. However, (qj,0) is not a saddle point 

or a minimum point with respect to the strains and, therefore, the deriva- 
tive at this point with respect to the strains does not vanish: 



CF-2 back 4926 


KA.CA oaf 4408 


11 


Equations (22) and (23) therefore become 


<P <P(<§,0) + (ag^q^ 0j0 8q J 8q k + 2^ ^a^ & r ^)o,0 


I 

d>»>P 


3>k 


fefe)o,o 69 ^ + • • • 


a,P 

r, s 


and 


<P = <p(q.f. 


’ o) + s tel, •<* + (wfe)*, 

a., P j , k 


8q,5q - 

,0 J K 


*Z 

a,p 

r^s 


a 2 ®. 


^ e cLp^ e rS/ ^ r8 


e~ D e-* + y| f v r — S( 1 a 6 „ 0 + . 


W 5 a8° q Jt,0 J 




Substituting equations (27) and (28) into the denominator 
tor of equation (20 ) , respectively, gives 


r 

s 



g(s) 


j f^(q) exp 

’ 1 V 5 , e / & 2 9 \ 

KE ^ “P \5eop5qJ 

j^,p 

/ f Q (q) exp 

Jv 

1 \ ' 8n c f. d 2 <? \ 

3>a-,P 


ap 6 ^ + 


(27) 


. . (28) 


and numera- 


dc 

0 


dv 

0 


( 29 ) 
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where r g is the jump frequency in the strained system and where the 
functions g( e ) , f-^q), and f Q (q) are defined by 


g(e) = exp 


*, 0 £ “ p ■ 

a, 8 a, 8 


s„ Q e, 




f a*; \ 

w s ap° 6 r sy o, o 


a, 8 

r,& 


\° c a8° s rs/t,o 


(30) 


and 


f*(a) - exp 


and 


f 0 (q) = exp 


“ M 5<1 J 5q k 

j> k 

“ kT ^<5'°) ‘ gkjf^^ (^5^) 0 , 0 6<1 J 6<lk 

Tk 


(31) 


(32) 


The functions f-j-(q) and fo(q) are the Boltzmann factors of the 

potential energy expanded about the saddle point of the activated state 
and about the minimum of the normal crystal. Therefore, in the unstrained 
case, equation (20) may be written as 


r u \2* 


i 1 / 2 


j f*(q) &<* 

f f o^) dv 

Jv 


(33) 


Therefore, multiplying and dividing equation (29) by T u , the jump 
frequency in the unstrained crystal, gives 


exp 


r o - r u g(«) 


kT ^ ^ 8q -j £ a8 (ss^Sqjjt^o 

J>a,P 


( exp - W ^ 8q ^ 

\ L vTs 


KT ^ -J 5 ae Rsjj^JI 0 0 

3,a,t 


(34) 
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where the statistical mechanical averages are given by 



(36) 

If the exponentials in equation (34) are expanded and only the first 
two terms in each expansion are retained. 
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However, in both the activated and normal states, the average dis- 
placement of. an atom from its mean position is zero; that is, 

(Sqj) a = ^)v =» 0 

and equation (37) reduces to 

r = r u g(e) (38) 

If the quadratic terms in the strain in equation (30) are neglected 
relative to the linear terms, substituting equation (30) into equation 
(34) gives 



Equation (39) shows that the jump frequency has a simple exponential 
dependence on the strains and that this dependence Is controlled by the 
derivatives of the potential energy with respect to the strains evaluated 
at the saddle point of the activated state. 

Equation (39) gives the general relation between the jump frequency 
and the strain that will be used in this report. If the strain matrix 
and the Interatomic forces are known, the effect of the strain on the 
jump frequency can be computed. 

To illustrate the application of equation (39), three special cases 
will be considered: 


(1) Uniform compression or expansion In which 

£ xx = 6 yy “ s zz “ 6 

(all other strains = 0) 

(2) Simple shear In which 

s xy = 6 yx *T e s 

(all other strains = 0) 


(40) 


(41) 
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(3) Simple elastic tension or compression in the x- direction, in 

•which 


xx 


®yy “ 6 zz “ -^®T 


where p. is Poisson’s ratio 

For these three cases, equation (39) gives: 
For uniform compression or expansion. 


r ° ‘ r “ « [- i (4,o 


For simple shear. 


T s = T u exp 




For simple elastic tension or compression in the x- direction. 


r s = r u exp 


®L \ , 2 ^ s T ( 5<p \ 


(42) 


(43) 


(44) 


(45) 


Since strong repulsive forces come into play as the diffusing atom 
moves to the activated position, the derivatives in equation (39) and 
equations (43) to (45) are all negative. Thus, negative strains, cor- 
responding to a compression of the lattice, decrease the jump frequency 
and positive strains, corresponding to an expansion of the lattice, result 
in an increase in the jump frequency. This conclusion is in accord with 
what is expected on a simple physical basis . 


The preceeding equations are in agreement with the results of refer- 
ence 9. On the basis of a simple model, which takes into account only 
nearest neighbor interaction, the effect of internal strains resulting 
from impurity atoms on the diffusion coefficient was computed (ref. 9), 
and it was found that the diffusion coefficient is an exponential func- 
tion of the strain. 


DEPENDENCE OF VACANCY CONCENTRATION ON STRAIN 


It was pointed out in the section BASIC EQUATIONS OF DIFFUSION THEORY 
that the quantity n appearing in the generalized Fick’s laws (eqs. (l7) 
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and (19)) lias a different interpretation for different mechanisms of dif- 
fusion. For dilute interstitial diffusion, n =* 1 provided there are no 
sources of interstitials, and Vn = 0 whether or not the system is 
strained. For diffusion by a vacancy mechanism, however, n is the atomic 
fraction of vacancies in the crystal given by 



(46) 


where 

n v vacancy concentration 

Nrp total number of lattice sites per cc. 

It is therefore necessary to investigate the variation of n y with strain. 

The concentration of vacancies in a crystal at equilibrium is given 
by 

n v " N T exl> (• AG v/ k3 ') (47) 


where 

n y vacancy concentration in the absence of strain 

AG V Gibb's free energy of formation of a vacancy 

In a strained system, the free energy of formation of a vacancy may 
be different from that in an unstrained system, so that the vacancy con- 
centration depends on the strain. 

If the strains are constant in time, the vacancy concentration is 
also constant in time. However, during plastic deformation, vacancies 
are produced at a rate that depends on the strain rate (ref. 10) . The 
general equation for the vacancy concentration is then 

n® = ^ + W T exp (-AG 8 /kT) (48) 

where 

n excess average concentration of vacancies arising from the plastic 
X strain 

AG 8 Gibb's free energy of formation of a vacancy In the strained system 
v 


4926 
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The analysis in this section consists of two parts: (l) the effect 

of static strains on the free energy of vacancy formation, and. (2) the 
effect of plastic flow on n x . 

The Helmholtz free energy A of a perfect crystal is given by sta- 
tistical mechanics as 

( 

I exp ( 


3N 


-to/ KE ) 77 a 2j 

3=1 


(49) 


exp (-A/kT) = 


for a system of N particles, where t^q is the energy of the system, 
expressed as a function of the coordinates q.j and momenta p ^ . For a 

crystal containing a single vacancy, the Helmholtz free energy is given 
by 

f [ 3N 

exp C-Ay/kT) = I ... I exp (-+ v /kT) TT dpj dq.j (50) 

J J J - 1 

where +0 and tjr v are the energies of the system in the perfect and 

Imperfect crystals, respectively. Therefore, since AG = AA + P AV, the 
Gibb's free energy of vacancy formation is given by 


exp (-AG v /kT) = 


/-/ 

hi 


3N 


exp (-ty/^T) “|"j" dpj dq^ 


, 1=1 


3N 


exp (-P AV /kT) 


exp (-t 0 /dT) JJ dpj dq.j 
3=1 


(51) 


where 

P pressure 

AV v volume change upon formation of a vacancy 
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If the phase integral is separated into configuration and momentum 
integrals, assuming the classical statistics of a crystal vibrating with 
normal mode frequencies, equation (5l) becomes 


exp (-AG V /W?) = 


Jexp (-qv/kT) a q TT < v 0>j 

j exp (-qi 0 /KC) dq "P[" (v T )j 


exp (-P AV /kT) (52) 


where (v q) j is the vibration frequency of the j vibrational mode in 
the perfect crystal and (v v )j is the frequency of the vibrationaJL 

mode in the crystal containing a vacancy. The integrals in equation (52) 
are for the entire configuration space. 


In a strained crystal, the free energy of vacancy formation is given 
by an equation analogous to equation (52) 


exp (-AG®/KP) 


Jexp (-<p|/ kT) dq (v 0 )j 

J exp (-<Po/kT) dq "j"j" (v®)^ 


exp (-P AV^ v /kT) 


(53) 


where the index s refers to the strained crystal. 

An estimate of the effect of strain on the vibration frequencies can 
be made from Gruneisen's relation (ref. ll) 


d In v _ 
d In V “ 


(54) 




where 

V volume 


X positive constant 
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Integrating equation (54) for each of the vibrational modes as the 
crystal goes from the unstrained to the strained state yields . 


<" 5> 4 / ' AvW 

V 1+ ^7 


(vfj 
I 


a ■ (> *r 


(55) 


It is evident that Gruneisen’ s relation leads to an equality of the 
products of frequency ratios in the strained and unstrained systems, so 
that combining equations (52) and (53) gives' 


exp 


- ^ (AG® - AG V ) 


) = 


I- 


I 


exp (-<p v /kT) dq I exp (-cp Q /kT) dq 


/ 


■ 


exp 


exp (-q>Q/kT) dq / exp (-qy/kT) dq 


[-(P AV^ - P AV v )/kr] 
(56) 


The volume difference (AV^. - AV y ) is given by 

(AV^ - AV y ) = (V® - V s ) - (V T - V) 
= (V^ - V v ) - (V B - V) 


where 


V 

volume 

V s 

volume 

V v 

volume 

V s 

V 

volume 


of perfect unstrained crystal 
of perfect unstrained crystal 

of unstrained crystal containing a vacancy 

of strained crystal containing a vacancy 
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For small strains, the volume change resulting from the strain should 
be relatively insensitive to whether or not the crystal contains a vacancy. 

The volume differences (V^. - V y ) and (V 0 - V) are therefore very nearly 
equal and it is sufficiently accurate to take the difference (AV^. - AV v ) 

equal to zero, so that the pressure-volume terms can be dropped from 
equation (56) . 

If the potential energy functions <p® and <p ® are now expanded in 

terms of the strains and the atomic displacements about the set of points 
corresponding to zero strains and mean atomic positions in the imperfect 
and perfect crystal, respectively, results analogous to equations (22) 
and (23) are obtained. Continuing the analysis leading to equation (39) 
shows that, to the first order, the free energy of vacancy formation 
depends upon the derivatives of the potential energy with respect to the 
strains, evaluated at the points in configuration space corresponding to 
the atoms in their mean positions in the perfect and imperfect crystals. 
However, the states corresponding to the perfect and imperfect unstrained 
crystals when 1 the atoms are in their mean positions are both equilib- 
rium states in which the systems are in potential wells with respect to 
both atomic coordinates and strains. Thus, the first derivatives are all 
zero and the free energy of formation of a vacancy is independent of the 
strains to a first-order approximation: 

AG® “ AG 

V y 

The terms quadratic in the strains are, of course, not zero. 

The analysis shows, therefore, that the jump frequency is more sensi- 
tive to strains than the energy of vacancy formation. In any diffusion 
process,, the effect of strain is felt most strongly through the jump fre- 
quency, and for static strains this is usually the only factor that must 
be considered. Thus, to the first order in the strains, equation -(48) 
may be written in the simple form 

n® = n x + Hqi exp (-AG v /kT) 

The only problem remaining is to calculate the effect of plastic flow on 


Theory and experiment both indicate (refs. 2, 10, and 12) that for 
simple types of deformation, the number of vacancies produced by plastic 
strain is proportional to the strain rate . In this report only simple 
plastic deformations that can be described by a single strain parameter 
are considered, Including tension, compression, and shear. 
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Accordingly, for the production rate of vacancies during plastic flow it 
is assumed that 

n + = %6 (57) 

where 

s strain rate 
Kj a constant 

Since the vacancy concentration during plastic flow is greater than 
the thermodynamic equilibrium value, the excess vacancies will tend to 
precipitate out of the crystal matrix. It has been shown that vacancy 
condensation is most probably a heterogeneous nucleation process (refs. 

13 and 14) in which the vacancies precipitate at imperfections in the 
crystal such as grain boundaries, foreign inclusions, and voids. For such 
a mechanism it is reasonable to postulate that vacancies are removed from 
the lattice by a first-order process. That is, the rate of destruction 
is given by 

n_ = Kgr^ (58) 

so that the differential equation governing the vacancy concentration is 

dn x , _ 

“dt = K l 6 " K 2 n x C 59 ) 

where 
t time 
Kg a constant 

On a microscopic scale, it is obvious that the excess concentration 
of vacancies may vary considerably from one point in the crystal to 
another, depending upon the distribution of the imperfections in the lat- 
tice, which act as vacancy sources and sinks . However, in macroscopic 
diffusion experiments, interest lies in the over-all average vacancy con- 
centration. Accordingly, a bar is placed over n^. to indicate the space 
average of the vacancy concentration. 

Integration of equation (59) gives 
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and equation (48) becomes 

4 = [l - exp (-K 2 t)] + % exp (-AG y /kT) (6l) 

The mode of production of vacancies by moving dislocations is dis- 
cussed in some detail in references 10 and 12. As pointed out in these 
references, a number of possible mechanisms exist by •which moving dislo- 
cations can generate vacancies. At present it does not seem possible to 
perform an accurate calculation of the rate of production of vacancies 
and therefore Kp must be treated as an empirical parameter to be deter- 
mined by experiment. However, it can be noted that if the generation 
mechanism involves dislocation climb, Kp will be proportional to the 

coefficient of self-diffusion as well as to the dislocation density. If 
a purely geometric mechanism is involved, Kj_ will depend only upon the 

dislocation density. Thus, for production by a geometric mechanism, 

and for production by a thermally activated mechanism, 

% = c 2 exp (-E/H) 


where 

c lt c 2 temperature independent constants 

E energy of activation for self- diffusion 

The rate of destruction of vacancies is determined by K 2 . Refer- 
ences 13, 15, and 16 are concerned with calculating the rate of absorption 
of vacancies by various types of vacancy sinks, and the remaining discus- 
sion of this section is largely an extension and development of this work. 

The constant K 2 depends on the nature of the vacancy sinks in the 

material. Obviously, different kinds of vacancy sinks are possible. 

Thus, a foreign Inclusion may collect vacancies and give rise to a void 
that acts as a spherical or ellipsoidal sink, and large grain boundaries 
may act aB two-dimensional sinks. Dislocations may act as sinks for 
vacancies in two ways: If the energy of interaction between a vacancy 

and a dislocation line is great enough, a vacancy becomes immobile and is 
essentially removed from the diffusion region when it gets closer to the 
dislocation line than some specified distance R Q . The dislocation line 

then giveB a rise to a cylindrical sink of radius R 0 « This type of sink 
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is proposed in reference 16. If the energy of interaction between a 
dislocation and a vacancy is large only at certain discrete points along 
the dislocation line, such as dislocation jogs, then K 2 depends on the 

frequency of co lli sion between a vacancy and a jog. This type of sink 
is postulated in reference 13 in studies of precipitation of vacancies 
during diffusion. Thus, four types of idealized vacancy sinks are con- 
sidered: spherical, platelike, cylindrical, and discrete-point sinks. 


Spherical Vacancy Sinks 


It is assumed that a spherical sink of radius R Q is located at the 

center of a spherical region of radius R. During plastic flow, the pro- 
duction rate of vacancies is K-j_£, so that the spherically symmetric 

steady- state diffusion equation is 


J v d 
„2 dr 



0 


(62) 


where 

D v diffusion coefficient for vacancies 

r radial distance from the center 

o steady- state conditions 

For the purpose of this calculation, the effect of the deformation 
on the coordinates in equation (62) Is neglected and It Is also assumed 
that D y is independent of position. 

The boundary conditions chosen for the solution of equation (62) are 

n£(R 0 ,t) = 0 


= 0 

r=R 

Equations (63) imply that the sink is a perfect absorber of excess 
vacancies and that the sinks are uniformly distributed through the crystal 
In such a way that the distance between them is 2R. 
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The solution of equation (62) with the boundary conditions given by 
equations (63) is 


o £lf /r3 r3 Bo r2\ 
n x “ SDy. ^R q “ r + 2 " 2 J 


(64) 


The average concentration n^. throughout the spherical volume of 
radius R is given by 




f 

Jro 


n£r 2 dr 


(R 3 - R 3 ) 

Performing the integration in equation (65) gives 


(65) 


-o K l e 

*x = “IT 


r 3 R 3 (R 2 - R 2 ) R 2 x (R 5 - Rq) 

+ 


3R, 


o 2(R 3 - R 3 ) 6 10 (R 3 - R 3 ) 


( 66 ) 


If R 0 < < Rj so that the distance between sinks is much larger than 
their radius, then, at steady state. 


-o _ fdf rJ 
n x " 3D V R 0 


But from equation (59), the steady-state condition gives 

Ki'e = K 2 n° 

which, combined with equation (67), gives 


(67) 


( 68 ) 


*2 = 


3 ®o -^v 


( 69 ) 


or, if there are f sinks per unit volume. 

Kg — 4 if R 0 fD v 


(70) 
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Platelike Sinks 

For platelike sinks of thickness Lq and a distance 2L apart, the 
boundary- value problem analogous to equations (62) and (63) is 

= 2 o 
d n* 


D y — g + K l 6 “ 0 


dx" 


^(^o) - o y 


(<) 

\ 


■ 0 


(71) 


J 


where x is the perpendicular distance from the plate. The solution of 
equation (71) is 


n S “ TC [O - “o) - I •- L o>] 


(72) 


and if L Q < < L, the average concentration is 


-o K 1® L 
n X " 3Dy 


(73) 


2D, 


K 2 = 


(74) 


or, if the area of the platelike sinks is approximately L 2 (as would be 
the case in a material of uniform grain size), then. 


K 2 * 2KyfL 


(75) 
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Cylindrical Sinks 

For cylindrical sinks* the boundary- value problem becomes 


55 ( r] V^S) +K i* 

nj(r 0 ) 


0 

0 


■ ° 


■where 

r Q radius of cylindrical sink 
2rj distance between sinks 

The solution of equation (76) is 

,2 


„o . K l e 




and for r < < r, , 
o l 7 


-o K l® 2 , r l 

n v = vrT r i ln ZT- 

x 2D 1 r 


and 


K 2 “ 


2D. 


v 


r? ln — 
- 1 - r_ 


(76) 


(77) 


(78) 


(79) 


If the sinks are dislocations of length l , the volume per sink is 
*r^7 and equation (79) can be expressed as 


*2 


2jtD y fZ 


(80) 
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But, fl is the number of dislocation lines per square centimeter, that 
is, the dislocation density Bj), so that 


k 2 = 


2*Dy%) 

i ri 

In — 
r o 


(81) 


Discrete Point Sinks 

If the vacancies are destroyed by an atomic collision process, as 
would be the case if only certain points, such as dislocation jogs, are 
effective in trapping vacancies, the rate of vacancy destruction is pro- 
portional to the collision frequency between vacancies and jogs. 

The jump frequency of a vacancy is roughly D v /X^, and if Cj is 

the dislocation jog concentration, the collision rate between excess 
vacancies and jogs is 






(82) 


where Cj/% is the probability that when a vacancy jumps, it runs into 

a jog. If It Is assumed that all collisions result in a destruction of 
the excess vacancy. 


*2 = 


D v C j 

X 2 Ntjt 


(83) 


Inspection of the various expressions for K 2 shows that Kg is 
proportional to the diffusion coefficient for vacancies and to ,the con- 
centration of vacancy sinkB, and depends on the geometric c har acter of 
the sink. 


In addition to the mechanisms involving migration to sinks, vacancies 
may disappear by combining with interstitials. Since interstitials are 
much more mobile than vacancies (ref. 16, ch. 5), only the diffusion rate 
of the Interstitials must be considered in the recombination process. 

The jump frequency of an interstitial is D^/X 2 where is the diffu- 

sion coefficient for the motion of interstitial atoms-. If the concentra- 
tion of interstitials is n^. Kg is given by 


K 2 = 


¥i 

x 2 % 


(83b) 
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Thus, if the vacancies are destroyed by recombination with Inter- 
stitials , Kg is proportional to the diffusion coefficient for inter- 
stitials rather than to the diffusion coefficient for vacancies. 

In an actual crystal, several of the production and annealing mecha- 
nisms may be operating simultaneously, in which case the constants 

and Kg are given by sums of the special cases described previously. 


DIFFUSION COEFFICIENT AS FUNCTION OF STATIC AND DYNAMIC STRAINS 

In the previous sections, the theoretical framework for an analysis 
of the effect of strain on diffusion has been developed. In this sec- 
tion, the previous results are combined to give the final, functional, 
dependence of the diffusion coefficient on the strains. 

From equation (19) it is evident that for an isotropic solid under 
homogeneous strain the diffusion coefficient for the flow of species N 
is given by 


D s = a*| n s r 0 


(84) 


where 

X B lattice parameter 

a constant that Is determined by the crystal structure 
The subscript s refers to the strained crystal. 

Four special, cases of the general equation (84) are considered: 

(1) Diffusion by an interstitial mechanism in the presence of static 
strain (strain rate, zero) 

(2) Diffusion by a vacancy mechanism in the presence of static strains 

(3) Diffusion by an interstitial mechanism during plastic deformation 
(strain rate not zero) 

(4) Diffusion by a vacancy mechanism during plastic deformation. 

The four diffusion coefficients corresponding to these four cases will 

be labeled D<i>, dM, D (i), and D^y), respectively. 

6 6 s e 
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For case 1, n g = 1 as discussed in the section DEPENDENCE OF VACANCY 
CONCENTRATION ON STRAIN and T s is given by one of the equations (39), or 

(43) to (45) . For the purposes of this discussion it is assumed that the 
deformation is a uniform lattice contraction, so that equation (43) gives 
the dependence of the .Jump frequency on strain. The extension to other 
types of strain is obvious. Using equation (43) and the fact that the 
strained lattice parameter is related to the unstrained lattice parameter 
by the factor (l + s) for case 1 yields 


D^g) = a\ 2 (1 + e) 2 T u exp 


5s /ckg\ 

hr 


(85) 


For case 2, the vacancy concentration is governed by equation (47), 
and to a first-order approximation is unaffected by strain. Since 
n s = n v /Ntji, and for the homogeneous deformation case being considered r g 

is given by equation (43), the diffusion coefficient for case 2 is given 

by 


= <xX 2 =? F (l + e) 2 esq? 

G Iliji 11 



( 86 ) 


Since plastic flow does not affect the number of interstitial sites, 

n = 1 even In case 3, and a result identical to that of equation (85) 
s 

Is obtained 


D W - ax 2 r u [1 + e(t)] 2 exp (87) 

A significant difference between equations (87) and (85) Is that in 
equation (87), the strain and, therefore, the diffusion coefficient, is 
a function of time. Also, in equation (85), the strain that must be con- 
sidered is the actual lattice strain and not the observed macroscopic 
strain. For plastic deformations, these two strains are not, in general, 
equal. 


For case 4, the vacancy concentration is given by equation (6l). 


Thus, again using equation 
4 gives 


(43) and the fact 


D<y) 

e 


aX 2 r. 


n 


[l+e(t)] 2 exp - gii. jl- exp C-Kat)]j 


that 

n s : 

( n v 

K-j_e 

NT 

k 2 % 




for case 


( 88 ) 
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For interstitial diffusion, the diffusion coefficient in the un- 
strained case is 

dC 1 ) = a\ z r n (89) 

and for diffusion by a vacancy mec h a ni sm in the unstrained case 


Therefore, equations (85) to (88) may be written as 



(93) 


These equations show how the effects of various types of strains on 
the diffusion coefficient can be taken into account for different diffu- 
sion mechanisms. The extension of this analysis to other diffusion mecha- 
nisms is completely analogous to the present development and will not be 
given here. 

It should be noted that in the. expression for D^, tlie strain rate 
in equation (93) is the macroscopic strain rate, although the strain 
e(t) is the microscopic lattice strain. 

COMPARISON OF THEORY WITH EMPIRICAL DATA 

The theory presented in this report makes a number of predictions 
that can be checked by existing experimental data. In this section, an 
analysis of the validity of the theory is made by comparing the theoreti- 
cal results with diffusion data. The equations developed thus far are 


Q7.P,' 
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not always in the form most suitable for comparison -with experiment. 
Whenever necessary, therefore, the equations will be transformed into a 
convenient form. 

The only literature data available for testing equation (9l) are 
concerned with the effect of hydrostatic pressure on the diffusion 
coefficient. 

The experimental data usu all y give the diffusion coefficient as a 
function of pressure, so that from the pres sure- volume relation of the 
material, the data can be obtained as a function of strain. Since 
pressure-volume data are generally given in terms of AY/V q , it is con- 
venient to express s as a volume strain (V D is the volume of zero 
pressure) . The strain s is 


X s ~ X AX 
6 X X 

so that, in terms of volume strain. 


_ 1 AY 

= 3 V 0 


(94) 


(95) 


for small strains. Therefore, using equations (94) and (95) and also 
the fact that 



changes equation (9l) to 



Therefore, it is evident that a plot of In D(g) (l + AV/y o )“^/3 
against AV/V Q should be linear with a slope m given by 



and an intercept given by In . 
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Several investigators have obtained data on the variation of the 
diffusion coefficient with pressure that is suitable for testing equation 
(97). Reference 4 presents data for the self-diffusion coefficient as a 
function of pressure for sodium, phosphorus, and mercury up to pressures 
of 12,000, 4000, and 8000 atmospheres, respectively. The self- diffusion 
coefficient of liquid gallium up to pressures of 10,000 atmospheres is 
given in reference 17. The self- diffusion coefficient for single crystal 
zinc up to pressures of 10,000 atmospheres for diffusion in the directions 
parallel to and perpendicular to the c-axis is determined in reference 1. 

The electric conductivities of silver chloride and silver bromide 
have been measured as a function of pressure up to 300 atmospheres (ref. 

18) . Since the conductivity is proportional to the diffusion coefficient 
of the silver ion by the Nernst-Einstein relation, the data of reference 
18 are suitable for testing equation (97) . 

Plots of the variation of the quantity log d(v) (l + AV/V 0 ) - ^/^ 
against AV/v o for the self -diffusion of sodium, phosphorus, mercury, 
and gallium are shown in figure 1. The quantities log (l+ AV/V C )~ 2 /^ 

for single crystal zinc were plotted againBt the fractional change in 
lattice parameter A\f\, since this is a more natural unit for discussing 
diffusion in anisotropic crystals and the linear compressions perpendicu- 
lar and parallel to the c-axis are available. The zinc data are plotted 
in figure 2. 

Figure 3 gives log g plotted against AV/V Q for silver chloride 
and silver bromide where R is the resistivity. The volume change AV/V 0 
is small enough for the pressure range considered so that (l + A V/V 0 ) -2 / 3 
does not appreciably affect the results and can be ignored. 

Compressibility data (refs. 19 to 23) were used to obtain the appro- 
priate value of AV/V 0 for zinc, sodium, mercury, silver chloride, and 

silver bromide. For gallium, AV/V 0 was computed from the data of refer- 
ence 22 assuming that the form of AV/V 0 as a function is the same as 
that for mercury. The values of AV/V 0 for white phosphorus were com- 
puted from data from reference 24 assuming that the variation of the 
fractional volume change with pressure has the same form as that observed 
in reference 25 for black and red phosphorus. 

In all cases, the available compressibility data were extrapolated 
to the diffusion temperature. 
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The linearity of the plots presented in figures 1 to 3 shows that 
the form of equation (97) is valid for those systems investigated within 
the probable inaccuracies of the experiments and the calculations. 

The slopes of the plots given in figures 1 to 3 are related to the 
interatomic forces through equation (97). If the repulsive potential is 
steep for a certain metal, that is, the atoms are "hard", then for a given 
atom displacement arising from a lattice strain, the rate of change of 
the potential energy with strain is much greater than for a similarly 
strained crystal containing "soft" atoms. Thus, (ckp/SsO^Q should be 

large for hard atoms and small for atoms having slowly varying potential 
functions. The quantity -l/k (dcp/de)^ q = mT has been calculated from 

the slopes of the plots of (l + A V/V 0 )“ 2 /^ plotted against AV/V 

and the absolute temperatures of the available diffusion experiments. 

Table I summarizes these values of m and mT. 

The lowest values of mT are those for mercury and gallium. This 
is to be expected since they have relatively "soft" potential functions, 
and being in the liquid state, their atoms are highly mobile so that 
they can adjust to the motion of the diffusing atom to give the lowest 
possible values to the interatomic forces . The values of mT for the 
silver halides are among the highest, and this is entirely reasonable in 
view of the steeply rising repulsive potential a silver ion meets as it 
migrates from one stable position to the next. In zinc, mT is much 
larger for diffusion perpendicular to the c-axis than for diffusion 
parallel to the c-axis . This Is in accord with the fact that the nearest 
neighbor distances are closer in the perpendicular direction, so that 
when an atom migrates to the activated position, the change in the Inter- 
atomic forces is greater than for a corresponding process in a direction 
parallel to the c-axis. Of all the solids listed, sodium has the lowest 
value of mT. The interatomic potential varies relatively slowly for 
sodium; in fact, recent calculations (unpublished MCA data) show that 
the potential well is so broad that the pairwise potential is repulsive 
to distances as far out as 1.4 times the nearest neighbor distance in 
solid sodium. Thus the low value of mT for sodium is In agreement with 
its interpretation In terms of the interatomic forces. 

The fact that mT is so much smaller for the liquid metals than for 
any of the solids including sodium is indicative of the difference in the 
mechanism of diffusion in liquids and solids. In a liquid, the atoms are 
not constrained to remain at lattice positions so that diffusion occurs 
by a cooperative process Involving the migrating atom and its nearest 
neighbors . Thus , the change in the interatomic forces can be kept, to a 
minimum throughout the diffusion process and, consequently, mT would be 
very low. 
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Dependence of Diffusion Coefficient on Strain Rate 

In the absence of strain, the diffusion coefficient D is given by 
the usual expression 

D = ccX^v n exp (-AG/kT) (99) 

For diffusion by a vacancy mechanism, n is the atomic fraction of vacan- 
cies in the crystal at equilibrium. In a crystal undergoing plastic 
deformation, the diffusion coefficient is 

D^ = <xX 2 n v exp ( -AG /kT) (100) 

As usual, the subscript s refers to the strained system. 


If it could be postulated that at the stress levels at -which plastic 
flow occurs at diffusion temperatures, the effect of the b train on the 
quantities X, v, and AG are negligible compared with the effect of the 
strain rate on the vacancy concentration, then to a good approximation the 
( v) / ( v) 

ratio D g /D v ‘ 1 would be given by 



( 101 ) 


This postulate is a reasonable one since it is well known that moving 
dislocations produce large numbers of vacancies. AIbo, plastic flow occurs 
by the displacement of large blocks of material as a result of dislocation 
motion, so that the microscopic strains defining the relative atomic 
positions are much smaller than the macroscopic strains. In fact, it is 
highly probable that the microscopic strains are always below the elastic 
limit of the material. 

Recent measurements (unpublished MCA data) of the effect of dynamic 
plastic flow on the rate of diffusion of hydrogen through nickel substan- 
tiate this hypothesis. The diffusion coefficient has been found to be 
independent of the state of plastic strain for tensile strain rates of 

0.02 to 0.4 hour Since hydrogen diffuses through nickel by an inter- 
stitial mechanism, any effect of strain on the diffusion coefficient must 
manifest itself through the quantities X, v, and AG. The hydrogen dif- 
fusion experiments can therefore be interpreted as indicating that the 
effect of plastic deformation on quantities other than the vacancy con- 
centration is negligible. 


’926* 
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For the simple types of dynamic strains discussed in the section 
DEPENDENCE OP VACANCY CONCENTRATION ON STRAIN, n B /n is obtained by di- 
viding equation (61) by equation (47), so that 


D^> 



exp 



( 102 ) 


Equation (102) shows that D g increases with time up to an asymptotic 
steady state at which 


d(y) K-i s 

1 + Vv 


(103) 


The time required to reach the steady state depends upon the value of 
Kg. In order to obtain an estimate of Kg, it will be assumed that the 

predominant type of vacancy sink is a cylinder around a dislocation line, 
so that equation (81) is valid: 




2 * 



(104) 


Typical, values of D v at approximately 1000° C are in the range of 10“ ^ 

"9 7 

to 10 centimeter squared per second and a value of about 10 for 

2 « N D /ln (r^/r^ has been given in reference 26, so that Kg is in the 

range of 1 to 10 -2 . Consequently, the factor |jL - exp (-Kgt)J reaches a 

value of 0.9 in a time somewhere between 0.04 and 4 minutes. For the ex- 
ample stated, a diffusion experiment should last at least several hours if 

-Q -7 

D t = 10 , and at least several minutes if D T is as high as 10 . 


The only data available in the literature on the effect of plastic 
deformation on diffusion is for the case of iron over a restricted range 
of strain rates (refs. 2 and 3). Although the data are not extensive, 
and some doubt exists concerning the absolute magnitude of the effect 
(ref. 3), a linear relation between the diffusion coefficient and the 
strain rate seems to be valid. This is in agreement with the steady-state 
equation (103) . In view of the scarcity of the data in this field, not 
much more can be said concerning the agreement of the theory with experi- 
ment. Additional experimental work in this area is highly desirable. 



36 


MCA TN 4408 


Volume of Activation 

From a comparison of equations (97) , (98), and (100), the free energy 
of activation for diffusion in a system subjected to hydrostatic pressure 
varies with the volume strain according to 

AG B = AG - m~^ kl (105) 

’o 


since v s * v for small lattice strains . Applying the definition of the 
activation volume for small lattice strains gives 


But, 


AV g = -mk3? 


c 


S(av/v 0 )' 


Up 



(106) 


(107) 


where 0 is the compressibility, so that the activation volume is given 

by 


AV g = mpkT 


(108) 


Table II presents values of the activation volume calculated from 
equation (108) at atmospheric pressure for those systems for which data 
are available. Figure 4 is a comparison of the activation volume of 
self-diffusion in sodium calculated from equation (108) with the activa- 
tion volume calculated in reference 4 from experimental data. The agree- 
ment is seen to be good. 


CONCLUSIONS 

A theory was developed that relates diffusion rates to the state of 
strain of the material. Fick's laws of diffusion were generalized to 
include the strain. The generalized equations differ from the ordinary 
diffusion equations in that the flux of diffusing material is proportional 
to terms containing the strain gradient as well as to terms containing the 
concentration gradient. In addition, a molecular-kinetic theory was 
developed that relates the diffusion coefficient to strain in terms of 
the atomic properties of the system. 
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The effect of dynamic plastic deformation on the diffusion coeffi- 
cient was investigated by considering the rate of production of vacancies 
by moving dislocations and the rate of precipitations of vacancies at 
vacancy sinks. The resulting equation states that the diffusion coeffi- 
cient is a linear function of the strain rate. 

Several predictions that can. be checked by existing experimental data 
may be made from the theory: 

1. Por diffusion as a function of hydrostatic pressure, the diffusion 
coefficient is an exponential function of the volume strain. 

2. The rate of change of the diffusion coefficient with strain is 
related to the interatomic forces . The relations is explicit enough that 
the variation of the diffusion coefficient with pressure can be interpreted 
in terms of the interatomic potential energy functions of the material. 

3. Por diffusion under hydrostatic pressure, the activation volume 
can be calculated from the compressibility and the rate of change of the 
diffusion coefficient with volume strain. 

4. Dynamic plastic deformation increases the diffusion rate, the 
diffusion coefficient being linearly related to the strain rate at steady 
state . 


In every case for which data are available, these conclusions are in 
agreement with experiment. 

The general framework of the theory provides a basis for understanding 
the effect of strain on diffusion in terms of the molecular-kinetic prop- 
erties of the system and should provide a valuable tool for comparing 
diffusion rates for different states of strain, as well as for investigat- 
ing the mechanism of diffusion. 


Lewis Flight Propulsion Laboratory 

national Advisory Committee for Aeronautics 
Cleveland, Ohio, Sept. 2, 1958 
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TABLE I. - COMPARISON OF VALUES OF mT 
FOR VARIOUS METALS 


Metal 

Temperature , 

T, 

°K 

m 

mT 

Sodium 

363 

27.9 

10,120 

Zinc (1) 

580 

86.6 

50, 200 

Zinc (M) 

580 

34.7 

20,100 

Mercury (liquid) 

303 

6.5 

1,970 

Gallium (liquid) 

303 

6.5 

1,970 

Silver in silver chloride 

573 

91.9 

52,600 

Silver in silver bromide 

573 

128.0 

73,000 


TABLE II. - ACTIVATION VOLUMES FOR SELF-DIFFUSION OF 
VARIOUS METALS AT 1 ATMOSPHERE CALCULATED 
FROM EQUATION (108) 


Metal 

Temperature, 

T, 

°K 

Activation volume, 
cc/g-atm 



Calculated 

Experimental 

Sodium 

363 

12.3 

12.3 (ref. 4) 

Phosphorous 

(white) 

314 

71.7 

30.0 (ref. 4) 

Zinc (i) 

580 

3.0 

4.9 (ref. l) 

Zinc (II) 

580 

8.3 , 

16.9 (ref. 1) 

Mercury (liquid) 

303 

.62 

.57 (ref. 4) 

Gallium (liquid) 

303 

.62 

.55 (ref. 16) 

Silver, in silver 
chloride 

573 

10.3 


Silver in silver 
bromide 

573 

13.7 
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Volume change, AV/V 0 

(b) White phosphorous at 314° K. - 

- 2/3 

Figure 1. - Continued. Variation of log D e fl + y~\ plotted 

against volume change AV/V 0 for a elf-dif fusion of various 
elements . 
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Figure 2. - Concluded. Variation of log D s plotted against 

fractional change in lattice parameter AA/1 q for self -diffusion in zinc. 
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(b) Silver bromide. 

Figure 3. - Concluded. Variation of log l/R plotted against volume change £V/V 0 for mobility 
of silver at 573° K. 
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Figure 4. - Activation volume plotted against pressure for self-diffusion 
of sodium at 362° K. 
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